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Abstract 

We propose and investigate an exactly solvable model of non-equilibrium Lut- 
tinger liquid on a star graph, modeling a multi-terminal quantum wire junction. 
The boundary condition at the junction is fixed by an orthogonal matrix §, which 
describes the splitting of the electric current among the leads. The system is driven 
away from equilibrium by connecting the leads to heat baths at different tempera- 
tures and chemical potentials. The associated non-equilibrium steady state depends 
on § and is explicitly constructed. In this context we develop a non-equilibrium 
bosonization procedure and compute some basic correlation functions. Luttinger 
liquids with general anyon statistics are considered. The relative momentum distri- 
bution away from equilibrium turns out to be the convolution of equilibrium anyon 
distributions at different temperatures. Both the charge and heat transport are 
studied. The exact current-current correlation function is derived and the zero- 
frequency noise power is determined. 
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1 Introduction 



The universal features of a large class of one-dimensional quantum models, exhibiting 
gapless excitations with linear spectrum, are successfully described [HE] by the Tomonaga- 
Luttinger (TL) liquid theory [3]- [6]. This theory^ applies to various systems, including 
nanowire junctions and carbon nanotubes, which are available nowadays in experiment 
[10] -[13]. For this reason the study of non-equilibrium phenomena in the TL liquid phase 
attracts recently much attention [I3]-[18|. 

A typical non-equilibrium setup, considered in the literature, is the junction of two or 
more semi-infinite leads with electrons at different temperatures and/or chemical poten- 
tials. The junction is an interval of finite length L, where the electrons injected from the 
leads interact among themselves. This interaction drives the system away from equilib- 
rium. Differently from the equilibrium TL liquid on the line, the non-equilibrium model 
defined in this way, is not exactly solvable. Nevertheless, it is extensively studied [H]-[T5] 
by various methods, including linear response theory, bosonization combined with the 
non-equilibrium Keldish formalism and perturbation theory. 

One of the main goals of the present paper is to explore the possibility to construct 
and analyze an alternative exactly solvable model for a non-equilibrium TL junction. 
Since the universal features of such a system are expected to manifest themselves in the 
critical (scale invariant) limit, it is natural to shrink the domain of the non-equilibrium 
interaction to a point, taking L — > 0. For a complete description of the critical regime 
it is essential to take into account all point-like interactions, which ensure a unitary 
time evolution of the system. These interactions can be parametrized by a scattering 
matrix § localized in the junction point, as shown in the multi-terminal setup displayed 
in Fig|T} Each lead contains a TL liquid, which at infinity is in contact with a heat 




Figure 1: A junction with scattering matrix § and n semi- infinite leads, connected at 
infinity to thermal reservoirs with temperatures Pi and chemical potentials /z;. 



reservoir with (inverse) temperature and chemical potential /ij. Our first step below 
is to show that there exists a non-equilibrium steady state (NESS), which describes the 
TL configuration in FigjTJ This state is characterized by non-trivial time independent 

1 For some more recent reviews we refer to 0-0 • 
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electric and heat currents, flowing in the leads. The scattering matrix § is implemented 
by imposing specific boundary conditions at the junction. It turns out that the boundary 
conditions, which describe the splitting of the electric steady current in the junction, lead 
to an exactly solvable problem. In fact, we establish the operator solution in this case and 
investigate the relative non-equilibrium correlation functions in the NESS representation. 

The TL theory has been introduced originally [3]- [6] for describing fermion systems. 
It has been understood later on [19]- [21] that the fermion TL liquid is actually an element 
of a more general family of anyon TL liquids^J which obey Abelian braid statistics. In 
this paper we explore the general anyon TL liquid, obtaining the conventional fermionic 
and bosonic ones as a special cases. 

From the two-point anyon correlation functions we extract the NESS distribution 
of the TL anyon excitations. In momentum space this non-equilibrium distribution is 
a nested convolution of equilibrium distributions at different temperatures and chemical 
potentials. As expected, the convolution depends on the scattering matrix §, which drives 
the system away from equilibrium. We investigate also the NESS correlators of the electric 
and energy currents, describing in detail the charge and heat transport in the junction. 
The zero-frequency noise power is deduced from the two-point current-current correlation 
function, whose exact expression in terms of hypergeometric functions is established. We 
prove the breakdown of time reversal invariance as well. 

The paper has the following structure. In the next section we construct non-equilibrium 
chiral fields in a NESS on a star graph modeling the junction. We derive here the non- 
equilibrium Casimir energy and the heat current and compare the latter with the con- 
formal field theory result. In section 3 we develop a non-equilibrium finite temperature 
operator bosonization procedure. We also establish the operator solution, subject to the 
current splitting boundary condition at the junction. We show that this condition covers 
two different physical situations, corresponding to a junction with and without charge dis- 
sipation. The non-equilibrium correlation functions are investigated in section 4, where 
the anyon NESS distributions are derived. The charge and heat transport as well as the 
noise are also studied there. Section 5 provides a concise outlook of the paper and con- 
tains some general observations. The appendix collects some results about the asymptotic 
properties of the anyon NESS correlators. 

2 Non-equilibrium chiral fields on a star graph 

The fundamental building blocks of bosonization away from equilibrium are the free mass- 
less scalar field ip and its dual ip. The fields <p and tp propagate on a star graph T, which 
is shown in Fig. [2] and models the quantum wire junction. The edges Ei are half-lines 
and each point P in the bulk r \ V of T is uniquely determined by its coordinates (x,i), 
where x > is the distance to the vertex V and i = 1, ...,n labels the edge. Besides the 
massless Klein-Gordon equation, the fields <p and <p satisfy the duality relations 

d t <p(t,x,i) = -d x (p(t,x,i) , d x <p(t,x,i) = -d t p(t,x,i) . (2.1) 

2 Similar results hold [551 H3] for the Lieb-Lineger and Calogero-Sutherland models. 
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Figure 2: A star graph T with n edges modelling the junction of n quantum wires. 

The initial conditions are fixed by the equal-time canonical commutation relations 

[tp(t,x,i) , tp(t,y,j)]_ = [tp(t,x,i) , tp(t,y,j)]_ = 0, (2.2) 

[(d t <p)(t,x,i) , tp(t,y,j)]_ = [(d t tp)(t,x,i) , <p(t,y,j)]_ = -iS^x-y). (2.3) 

In order to determine the dynamics completely, one must impose some boundary condi- 
tions at the vertex x = 0. These conditions are conveniently formulated in terms of the 
combinations 

Wi,n{t - x) = (p(t, x, i) + (p(t, x, i) , <fi,L(t + x) = (f{t,x,i) - (f(t,x,i) , (2.4) 

which depend on t — x and t + x respectively and define right and left chiral fields (p^z 
on T. The most general scale invariant boundary conditions, generating a unitary time 
evolution of (p and (p, are parametrized by the orthogonal group 0(n) and read [24]- [27] 

n 

ViAO = E ' S G 0(n) . (2.5) 

These simple conditions capture the universal features of the system and S has a straight- 
forward physical interpretation: the vertex V of T represents a scale invariant point-like 
defect, § being the associated scattering matrix. 

2.1 The non-equilibrium steady state ^ ^ 

Our next step is to construct a steady state , which captures the evolution of the 
chiral fields (p^z on T, whose edges are attached at infinity to thermal reservoirs at inverse 
temperatures fa as shown in Fig. [TJ In the boson case we take all chemical potentials 
to be equaQ setting = /x 6 in all reservoirs. The system is away from equilibrium if S 
contains at least one non-trivial transmission coefficient among reservoirs with different 
temperature. The construction of , described below, follows the scheme developed 
in [28] and is based on scattering theory. It adapts to the case under consideration some 
modern ideas [29]- [33] about NESS. The framework is purely algebraic and generalizes 
the definition [31] of equilibrium Gibbs state over the algebra of canonical commutation 
relations (CCR). 

3 This choice will not prevent us to deal in the fermion case below with arbitrary jiij. 
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We start by observing that the massless Klein-Gordon equation and the relations (2.1 ) 
lead to the following representation 



<PiAO 



00 dA; 
00 dk 



(2.6) 
(2.7) 



/o vrv/2 

A\ being some distribution to be fixed below. Using that § is a real matrix, the boundary 



condition (2.5) implies the constraints 



i(*0 = ^2 s a(k)a j {-k) , a*(k) — ^*%(A;) a*(-k) 
3=1 3=1 



(2.f 



where 

S(k) = 9(-k)§ t + 9(k)§, (2.9) 

9 is the Heaviside step function and §* indicates the transpose of §. From the equal-time 
commutation relations (2.2, 2.3) one infers that the elements {cii(k), a*(k) : k e M, i = 



1, ...,n} generate the following deformation A of the algebra of CCR: 

[di(k) , a 3 -(p)]_ = [a*(fc) , a*(p)]_ = , 
[ai(fc) , a*(p)]_ = 2n[5(k-p)5 ij + S l3 (k)5(k + p)] . 



Moreover, (2.2, 2.3) imply that 



1*1 A A (fc) = 1 . 



(2.10) 
(2.11) 

(2.12) 



There exist a one-parameter family of tempered distributions, which solve this equation 
in R. A convenient representation of this family is given by [35] 



A x (k) 



dk 



9{k)ki 



(2.13) 



where A > is a free parameter with dimension of mass, having well-known infrared 
origin. 

The above structure is very general and equations ( 2.6|2.7 ) apply to any representation 
of the algebra A, which is a simplified version of the so called reflection-transmission 
(RT) algebra |36]-j38], describing factorized scattering in integrable models with point- 
like defects in one dimension. The Fock and the Gibbs state over A describe equilibrium 
physics and have been largely explored. We will investigate here the NESS f2g iM , which 
describes the physical situation shown in Fig. [TJ For this purpose we first observe that 
the sub-algebras Am and *4. ut, generated by the elements {cti(k), a*(k) : k < 0} and 
{cii(k), a*(k) : k > 0} respectively, parametrize the asymptotic incoming and outgoing 
fields. Accordingly, both A in and -4. ut are conventional CCR algebras; in fact the S(k + p) 



term in (2.11) vanishes if both momenta are negative or positive. It is worth stressing 
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that (2.8) relate Am with A ou t and that the whole RT algebra A can be generated via 
(2.8) either by Am, or by A on t- The main idea for constructing £lp }lib is based on this kind 
of asymptotic completeness property. Starting with an equilibrium state on A m , we will 
extend it by means of (2.8) to a non-equilibrium state on the whole algebra A. For this 
purpose we introduce the edge Hamiltonian and number operators 



dk 
2^ 



k\a*(k)a,i(k) , rii 



dk 
2^ 



a*(k)a,i(k) 



(2.14) 



which describe the asymptotic dynamics at t 
of Am- Defining 



-oo (i.e. before the interaction) in terms 



8=1 



we introduce the equilibrium Gibbs state over Am in the standard way [3 
polynomial V over Am we set 



(2.15) 
For any 



(2.16) 



where fcj < 0, pj < and Z = Tr (e - ^). All the expectation values (2.16) can be 
computed [M] by purely algebraic manipulations and can be expressed in terms of the 
two-point functions, which are written in terms of the familiar Bose distribution 



bi(k) 



-Pi[\k\-»b] 



-API-w,] 



(2.17) 



in the following way 

(a* j (p)a i (k)) Pt to = b i (k)6 ij 2'K6(k-p), (2.18) 

( ai (k)a*(p))^ b = [1 + biik^nSik-p) . (2.19) 

We stress that ( 2.18[2.19 ) hold on A- m , i.e. only for negative momenta. The common for all 
reservoirs chemical potential /i.^ < allows to avoid in ( 2.18[2.19 ) the infrared singularity 
at k = 0. We anticipate that //& has nothing to do with the fermion chemical potentials, 
appearing in the non-equilibrium bosonization procedure described in the next section, 
where the limit — > 0~ exist and will be performed. 

The next step is to extend (2.16 2.19) to the whole RT algebra A, namely to positive 
momenta. Employing (2.8) one finds 

(a*(p)ai(k))^ b = 2vr 



n 



9(—k)bi(k)8ij + 9(k) > § u b t (k) % 5(k - p) 



i=i 



9(-k)bi(k) §*,■ + 9(k)S ij bj(k) S(k + p) 



(2.20) 
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The expression for (ai(k)a*{p)) p ifl is obtained from (2.20) by the substitution 

b m^i+bm = 1 _ e _l m _, b] - (2.2i) 

The final step is to compute a generic correlation function. By means of the commu- 
tation relations (2.1012.11 ), this problem is reduced to the evaluation of correlators of the 



form (11™=! Ui m {h m ) 111=1 °i„(PiJW»> wm ch can be computed by iteration via 



M 



N 



m=l n=l 

M M N 

<W ^(aiiihja^ipjj)^ (J[ a im (k im ) a * Jn (PjJ)p,n b 



(2.22) 



m=l 



m=2 



n=l 



We would like to mention in conclusion that the use of the RT algebra A in the 
construction of represents only a convenient choice of coordinates, which has a 

simple physical interpretation in terms of scattering data and applies to a variety of 
systems [281 ES] with point-like defects. 



2.2 Energy density and energy transport in 

In order to illustrate the physical properties of it is instructive to investigate the 

non-equilibrium energy density and transport associated with the scalar field (p. The 
equations of motion imply the conservation 

d t 6 u (t, x, i) - d x 6 xt (t, x,i)=0, (2.23) 

of the energy-momentum tensor 

e tt (t,x,i) = ± -<p(%<p)]:(t,x,i), (2.24) 

6 xt (t,x,i) = i : [(d x ip)(d t <p) -<p(d x d t ip)} : (t,x,i), (2.25) 



where : • ■ ■ : denotes the normal product in the algebra A. The boundary condition (2.5) 
implies the Kirchhoff rule 

n 

^MW) = 0> (2-26) 

i=i 

which, combined with ( |2.23[ ), ensures energy conservation. 

The derivation of (6 tt {t,x,i))p^ b and (9 xt (t,x,i))i3 jfib is based on the expectation value 

(: ip(t 1 ,x 1 ,ii)ip(t 2 ,x 2 ,i 2 ) :)m = 

"°° dk r 

—A x (k)^5 ili2 b il (k) cos[/c(ti2 + x 12 ] + S ili2 b i2 (k) cos[k(t 12 - x 12 ] + 

n 

b h {k) §* iJ2 cos[k(t l2 + x 12 ] + Siii h A k ) § k cos[A;(t 12 - x 12 ]} , (2.27) 
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where x\i — x\ + x 2 - Plugging (2.27) in the definitions ( 2.24|2.25 ), one gets 

£*(P,Pb) = (8tt(t,x,i))p^ b = 

f°° dk v-^ 1 

§ii / — k cos(2kxMk) + V (<% + S?-) —2 Li 2 (e^ 6 ) , 

n 

Ti&to) = (e xt (t,x,i)) fi> , b = (^-SJ) 2^ Li 2 (e^) , 



i=i 



(2.28) 
(2.29) 



where Li 2 is the dilogarithm function. Eq. (2.28) determines the Casimir energy, whereas 
(2.29) describes the energy (heat) transport. Both are time- independent, thus confirming 
that we are dealing with a steady state. Notice also that the energy density is x-dependent, 
which reflects the breaking of translation invariance by the junction and is consistent with 
the conservation law (2.23). Since § is an orthogonal matrix, Ti((3,fib) obviously satisfies 
the Kirchhoff's rule (2.26). The energy density fib) can be written in the equivalent 
form 



— k cos(2kxMk) + — 2 Li 2 (e ft ^) - ]T (% - §?■) — 2 Li 2 (e^) , 

* j=i j 

(2.30) 

where the last term vanishes at equilibrium and describes therefore the non- equilibrium 
contribution to the Casimir energy. 

For a junction with n = 2 wires there are two one-parameter families 



cos 9 sin 9 
— sin 9 cos 9 



cos 9 sin # 
sin # — cos 9 



9 e [0, 2vr) 



of scattering matrices with det(§ ± ) = ±1. For both families one finds 

sin 2 



7i(/3, Hb) = -T 2 (f3,fi b ) 



2tt 



^ (e**) - ^Li 2 (e^) 



0i 



ft 



(2.31) 



(2.32) 



which deserves a comment. The heath transport has been investigated recently in the 
framework of conformal field theory (CFT) for generic central charge c in [3D]. Eq. (2.32) 
confirms the result of [30] for c = 1 and extends this result in two directions: imperfect 
junction with transmission probability sin 2 (6 l ) and //& 7^ 0. We observe in this respect that 
derivation of the dilogarithm terms in (2.32) is problematic in a CFT context, because 
both Pifib and 02 fib are non-trivial dimensionless parameters. 

In the limit fib — > 0~ the x-dependent integral in (2.28) can be evaluated explicitly 
and one finds 



W,0) 



7T 

12 

7T 

12 



n ^ ^ 



7T 



J'=l 



2)8? 



sinh ( 27r 



1 2 ' 



i=i 



(2.33) 
(2.34) 
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2.3 Chiral NESS correlators 



In the non-equilibrium bosonization procedure, developed below, we will need the corre- 
lation functions of the chiral fields ( 2.6|2.7 )) in the NESS £lp,n b - It is easily seen that all 
of them can be expressed in terms of the distribution 



dk 



71 



Aa(*) 



e ifc « + 



1 



1 — Q-P[k-Hb] 



The full A-dependence and the singularity at fib = of (2.35) are captured by 

1 f 



71 



UN 



hi 



x 



In 



71 



2i sinh — £ — ie 
1 f3 



o(nb 



(2.35) 



(2.36) 



where o(fib) stands for A-independent terms, which vanish in the limit //& — > 0~. It is 
convenient at this point to relate the (up to now free) infrared regularization parameter 
A to fib by means of 

(2.37) 



A 



The limit fib ~ > in (2.35) now exists and gives the distribution^] 

1 



w(£,p)= lim w(£, (3;\ = \fi b \,fi b ) = In 

/Lii,->0- 71 



2i sinh 



7T 







£-ie 



(2.38) 



which is the fundamental block of the NESS chiral correlation functions. In fact, for the 
two-point correlators one gets 

(<Pii,Lfa)'PizA&))p = S hi2 W (€l2,Ph) , (2.39) 

((PiuLfofoitA&h = ^(62,/3n)SU, (2.40) 

(<^i,i?(£iVi 2 ,£(6))/3 = S ili2 w(^2, A 2 ) , (2.41) 



(^i,-R(Ci)^i 2> fl(6))^ 



J«2 



(2.42) 



^§n^(ei2,/3,)§5 

where £12 = £1 — £2- As expected, the point-like interaction in the vertex of T induce a 
non-trivial left-right mixing described by Q2.40|2.41|). 



3 Bosonization away from equilibrium 

The possibility to express fermions in terms of bosons in 1+1 dimensional space-time has 
been discovered long ago by Jordan and Wigner [JT]. The bosonization technique in the 
Fock representation of the fields ip and tp has been applied for solving the Tomonaga- 
Luttinger (TL) model in [lj-[6j. The framework has been extended later (T9J H2] to the 
finite temperature Gibbs representation of ip and <p. Both the Fock and Gibbs represen- 
tations describe equilibrium physics. Our goal in what follows we will to apply the NESS 
representation, constructed in the previous section, for investigating the non-equilibrium 
TL liquid in the multi-terminal configuration shown in Fig. [TJ 

4 The ie prescription, adopted throughout the paper, indicates as usual the weak limit e — > + . 
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3.1 The Tomonaga-Luttinger model on T 

The bulk dynamics is governed by the TL Lagrangian density 

£ = iifi(d t - v F d x )^ + ty*(d t + v F d x )^ 2 - g+^Wi + - g-tyWi ~ , (3.1) 

where {ipa{t, x,i) : a = 1, 2} are complex fermion fields, v p > is the Fermi velocity and 
g± G K are the coupling constants. 

The bulk theory has an obvious Ul(1) <g> Ur(1) symmetry. In fact, the Lagrangian 
density (3.1) is left invariant by the two independent phase transformations 

Tp a ^e is <*ip a , C^e-^C, s a ER, a = 1,2. (3.2) 

implying the current conservation laws 

d t pz(t, x, i) - v F d x j z (t, x,i) = 0, (3.3) 

where the charge and current densities are given by 

/[^ViJ^^i), Z = L, f^JVi] Z = L > 

\[^ 2 ](t,x,z), Z = i?, JZ{ ' ,J \-[ftih](t,x,i), Z = R. 

(3.4) 

The currents j z (t,x,i) have simple physical meaning: j L (t,x,i) and jn(t,x,i) represent 
the particle excitations moving along the edge Ei towards and away of the vertex V 
respectively. Interpreting the vertex as a defect, which can be characterized by some 
scattering matrix, the currents jh and jr describe therefore the incoming and outgoing 
flows. 



3.2 The current splitting boundary condition 



It is well known [T]-[6] that the TL model (3.1) is exactly solvable on the line M. On the 
graph T the situation is more involved, because one should take into account the boundary 
conditions in the vertex V. The conditions 



n 

E 



U G U(n) 



(3.5) 



which work in the free case g- = g + = 0, do not lead [13] to exactly solvable problem after 
switching on the TL interactions. The fact that on R the quartic bulk interactions in (3.1 ) 



are solved exactly via bosonization suggest to try boundary conditions which, differently 



from (3.5), are formulated in terms of real boson fields. In this spirit and according our 



previous comments on the chiral currents (3.4), it is quite natural to consider 



j R (t,0,i) = J2 S ik3L(t,0,k), JeO(n) 



(3.6) 



k=l 
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which has been proposed and explored first in the two-terminal case in [H] . An advantage 
of (3.6) is the direct interpretation in terms of gauge invariant physical observables, which 
represent the basic building blocks of algebraic quantum field theory (see e.g. |45J). In 
fact, (3.6) describes the splitting in the vertex V of the outgoing current jn(t,Q,i) along 
the edge Ei in incoming currents ji{t, 0, k) along the edges E^. For this reason we refer 
to J the as the current splitting matrix and show in the next subsection that J actually 
coincides with the boson scattering matrix S. 



3.3 Operator solution of the TL model on T 

Referring for the details to jl6] , we recall here the anyon operator solution of the TL model 
on a star graph 17 The solution provides a unified description of all anyon Luttinger liquids 



and is expressed in terms of the chiral fields ( 2.6|2.7 ) and the parameters a, r G M. and 
the sound velocity v G K as follows: 



ip 2 (t,x,i) 



•q. ■ e i V^[ Tl Pi,R( vt - x )+ cr 'Pi,L( , "t+x)} . 



(3.7) 
(3.8) 



Here : ■ ■ ■ : denotes the normal product in the RT algebra A and r\i are some Klein 
factors, controlling the statistics of ip a . In this respect we impose the general anyon 
exchange relation 



i[)l(t,x 1 ,i 1 )i/; a (t,x 2 ,i 2 ) 



3 (-l) a i7rKe(a;i2) 



lf) a {t,X2,i2)i>a( t > x h'i>l) , X!^X 2 , (3.9) 



where e(x) is the sign function and k > is the so called statistical parameter which 
interpolates between bosons (k, - even integer) and fermions (k - odd integer). A simple 
realization of the Klein factors is 



1 



2ti 



: c 



Ti(7i+7*) 



where {7^, 7* 



1, ...,n} generate the auxiliary algebra 

fo > 7,] = [7* , 7*] = , h , 7*] = 



(3.10) 



(3.11) 



with 6y = —1 for i < j, en = and e^- = 1 for i > j. 

In order to fix the solution ( 3.7|3.8 ) completely, one should determine the parameters 
a, t and v in terms of coupling constants g± and the statistical parameter k. Using a 
standard short distance expansion and ( 3.7|3.8 ) one gets the charge and current densities^] 



P±(t,x,i) = (: Vi^i : ± : ^2^2 -){t,x,i) 



-1 



5 Without loss of generality we assume in what follows r > and r 7^ ±a 



[(d<p itR )(vt -x)± (d(pi, L )(vt + x)] , 

(3.12) 
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j±(t,x,i) 



2y/7TV F (± 

where for convenience the variables 



[(d<Pi,R)(vt - x) =p (d<f ijL )(vt + x)) 



C±=r±a. 



(3.13) 



(3.14) 



have been introduced. The normalization of (3.12) is fixed [46J by the Ward identities 
associated with the electric charge Q + and the helicity Q_ defined by 



n roc 

Q± = ^ / dxp±(t,x,2) . 

i=l - 70 



The normalization of (3.13) in turn is determined by the conservation law 

d t p±(t,x,i) - v F d x j±(t,x,i) = . 

Plugging ( 3.7|3.8|3.12 ) in the quantum equations of motion 

i[d t + (-l) a v F d x ]^ a (t,x,i) = 
%+ ■ p+(t,x,i)ilj a : (t,x,i) - (-l) a g- : p-(t,x,i)i/; a : (t,x,i)\ , 



(3.15) 



(3.16) 



one finds 



(3.17) 

(3.18) 
(3.19) 

(3.20) 

Eqs. ( 3.20|3.18|3.19 ) provide a system for determining v and C± ( or equivalently a and 
r) in terms of K and g±. The solution is 



7T 

2 

Vir/t H — g- . 



Moreover, the exchange relation (3.9) implies 



(I 



/ , o \ ±1/2 

/ TTKVp + 2g + \ 



K ~ — 

\tik,vf + 2g_ 
(ttkvf + 2g_)(irnvF + 2g + ) 



■kk, 



(3.21) 
(3.22) 



where the positive roots are taken in the right hand side. The relations (3.21) and (3.22) 
represent the anyon generalization [16] of the well known result for canonical fermions 
k = 1, where an alternative and frequently used notation [7j is 



92 = 2(g+ - g-) , #4 = 2(g + + 



K = C = C 2 



(3.23) 
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Considering the general anyon solution (3.21, 3.22), we assume in what follows that 
the parameters {k, g±} belong to the domain 

V = {k > 0, 2g± > -ttkv f } , (3.24) 

which ensures that er, r and v are real and finite. 



and 8. Expressing the chiral currents jz in terms of the chiral fields 



Let us discuss finally the current splitting boundary condition (|3.6|) and establish the 
relation between 
(fi,z, one finds 

1 

r 
i 

2' 



3R(t,x,i) 



K ./ 



C+j+)(t,x,i) 



v 

TTVp 
V 



d<Pi,R(vt - x) 

d<fi,L{vt + x) 



(3.25) 
(3.26) 



which, according to (2.5) satisfy the current splitting boundary condition (3.6), provided 
that 

J = 8 G 0{n) . (3.27) 



The symmetry content of the TL junction is strongly influenced by (3.27). The point 
is that in the presence of a defect the continuity equation (3.16) alone is not enough to 
ensure the electric charge conservation. A direct computation shows indeed that 



d t Q, 




(3.28) 



- ^2§ik ) (d<p kiL )(vt) . 
i=i J 

The independence of (p^L implies that the electric charge Q + is conserved if and only if 

n 

J2^ k = l, Vfc = l,...,n, (3.29) 



which, as expected, is equivalent to the Kirchhoff rule 

■n 







(3.30) 



Since 8 G 0(n), one infers from (3.29) that 8* satisfies (3.29) as well. Therefore, the 
electric charge Q + is conserved for those 8, whose entries along each column (line) sum 
up to 1. In geometric terms, these scattering matrices belong to the stability subgroup 
O v C 0(n) of the n- vector v = (1, 1, 1). An explicit parametrization of O v in terms of 
angular variables is given in |47j . 

Summarizing, the condition 8 G 0(n) guaranties the energy conservation in the TL- 
j unction. Concerning the electric charge Q + , one must distinguish two different regimes. 
Q + is conserved for 8 G O v . If instead 8 belongs to the complement O v = 0(n)\O v , there 
is an external incoming or outgoing charge flow in the junction and Q + is not conserved. 
The possibility to describe such imperfect junctions is a remarkable feature of the current 
splitting boundary condition (3.6). The physical details about the charge transport in 
the junction are discussed in section 4.2 below. 
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3.4 NESS representation and chemical potentials 



The crucial property of the operator solution ( 3.7|3.8|3.21|3.22[ ) is that it is universal, 
meaning that it applies for any representation of the chiral field algebra generated by 
(fi,z- The Fock and Gibbs representations have been largely studied and describe the 
equilibrium properties of the TL model on T. In order to explore the behavior of the 
Luttinger liquid away from equilibrium, we investigate below the operator solution in the 
NESS representation of the RT algebra A, constructed in section 2.2. 

The first step in this direction is the introduction of the fermion chemical potentials 



fj,i = k F - Vi 



(3-31) 



where kp defines the Fermi energy for k = 1 and Vi is the external voltage applied to 
the thermal reservoir in the edge Ei of Fig. [2] In what follows we keep kp fixed and 
vary eventually the gate voltages V^. As already mentioned, the boson chemical potential 
fib < has been introduced for avoiding some infrared singularities at the boson level 



and has nothing to do with /ij. In fact, in the chiral correlators (2.39 2.42) we already 
performed the limit //& — > 0~ 
a^, defined by 



In order to recover /Xj, following [19] we introduce the shift 



<Pi,L(Q ' ► (<W,l)(0 = PiAZ) 



Mi 



and, consistently with the boundary condition (2.5), 



ViAQ 1 — ► (<w,.r)(0 = - 



3=1 



(3.32) 



(3.33) 



The transformations ( 3.31|3.32 ) extend to an automorphism on the whole algebra 



generated by the chiral fields <p%z\ which is directly implemented in the operator solution 



(3.7|3.8|3.12|3.13). At this stage the TL correlation functions in the NESS are defined by 



(Oi[<fi lt z) ■ ■ ■ O k [(p ihtZ ])i3 >lt = (0 1 [a^ luZ ] ■ ■ ■ O k [a^i k ,z])p 



(3.34) 



In the rest of the paper we focus on the correlation functions (3.34), which capture the 



physical properties of the Luttinger liquid with the current splitting boundary condition 



(3.6) away from equilibrium. We will show in particular that (3.34) satisfy the Kubo- 



Martin-Schwinger (KMS) condition [34, 45J at equilibrium, which justifies the introduction 
of the chemical potentials by means of q3.32|3.330 . 



4 Non-equilibrium TL correlation functions 



4.1 Anyon correlators 



We derive here the two-point correlators of ip a (t,x,i) defined by (3.7, 3.8) in the NESS 



and discuss their properties. For this purpose we extend away from equilibrium the finite 
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temperature results of [19J. Using (2.22), for ip\ one finds 

(ipKh, xi, i)ipi(t 2 , x 2 , j)) p tfl = AijBij(t lt 2,xi t2 ;ij) x 



1 



sinh 



f^vt^ + x^) - i£ 



sinh 



ft 



(vtu + X12) - is 



n 



{vt l2 - x 12 ) - is 
1 



k=l 



ft 



sinh 



j^(vt 12 - x 12 ) - ie 



— sinh 
where 

ij 2tt \2i 

B i j(t 1 2 , x\ 2; /i) = e i ^ r ^ t,<1+x ' 1 ^ i_ ^* 2+:r2 ^^ +cr ^ 1 ' <1_x ' 1 ^fc= 1 n ikttk-(vt2~x 2 )Y.k=i s jkHk}}/(o-+T) 



(4.1) 



i7r 2 KeiJ /2 / l \ (^ 2 +r 2 )<5 1J + CT r(§ lJ +S* J ) 



The ^-correlator has the analogous form, 



(^2(ti,£i,i)^ 2 (*2,£2, = (4.1) with a^r. 



(4.2) 

■) _ 
(4.3) 

(4.4) 



The TL junction involves two types of ip\-ip 2 interactions. First, the Lagrangian (3.1) 
contains a ipi-ip 2 bulk coupling proportional to (g + — gJ). Second, the current splitting 
boundary condition (3.6) provides an additional boundary interaction described by the 
mixed left-right correlators (2.40, 2.41)). Consequently, the mixed ipi-ip 2 correlators are 
non-trivial and have the form 

(-0* (ti,xi, i)?p 2 (t 2 , x 2 , j))^ = A ij B ij (t 1:2 ,x lj2 ;n) x 



— sinh 



1 



r 2 S'. 



sinh 



f-(vt 12 - x 12 ) - is 



ft 



in, fcj 



sinh 



where 



A; 



n 

^=1 



ft: 



sinh 



ft 



.vfo - X12) - ie 



(4.5) 



e i7T 2 Ke ij /2 / -|_ \ <T 2 § lJ +T 2 S,l J +2(TT8 i] 



2tt V 2i 



(4.6) 



B i j(t 1 2, X\ 2\ /J-) = e i ^ r ^^ <1+X ' 1 ^ i_ ^ i2 ~ X2 ^^ lS3ft ' lfc ^ +Cr ^ l '' 1_:ri - ) ^*= 1 n ik^k-(vt2+x 2 )^}}/(a-+T) 

Finally, 



(ip2 , (ti,x 1 ,i)ipi(t2,x 2 ,j))p,v, = (4.5) with cr-B-r. 



(4.7) 
(4.8) 
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As expected, in the equilibrium limit f3i = f3 and //, = // for all i, the correlators 



(4.1)-(4.8) simplify and satisfy the KMS condition, which represents a non-trivial check 



both on the computation and on the shift (3.32, 3.33) introducing the chemical potentials. 



Let us consider for instance (4.1), which in this limit takes the form 



1 



sinh 



|(u*i2 + xi 2 ) - is 



1 



<rr§' 



sinh 



|(ftl2 - X12) - ie 



1 



a 2 &i 



sinh 



I (vti2 + X12) ~ ie 



sinh 



|(vti2 - £12) - ie 



Recalling that the KMS automorphism £> s acts on tp a as follows, 



[Qs^a] (t, X, l) 



1pa(t + s/v,X,l) 



one can check that the equilibrium correlator (4.9) satisfies the KMS condition 



(V>J(ti,a:i,i) [Qs+iM (t 2 ,x 2 ,j))p,p = ([Qs^Pi] (t 2 ,x 2 ,j)4>* 1 (t 1 ,x 1 ,i)) /3 ^ 
for all values of the statistical parameter k. 



(4.9) 

(4.10) 
(4.11) 



The critical scaling dimensions d{ can be extracted from (4.1 )-(4.8 ) in the limit j3{ — > 00 



and ^ — > 0. Because of the operator mixing, this is a subtle issue, which has been 
discussed in full detail in |32]. One gets, 

1 



di 



(a 2 + r 2 ) 



GTS; 



, n . 



(4.12) 



where Sj = ±1 are the eigenvalues of §. As already observed in the impact of the 
vertex interaction is captured by the term arsi, which preserves unitarity in the sense of 
conformal field theory because d, > 0. 

A remarkable special case is obtained by setting g + = g_ = g. In this case the bulk 
ipi~ip 2 coupling vanishes and one is left only with the boundary interaction induced by the 



current splitting boundary condition (3.6). From (3.22) one gets 

2g 



V F + 



7TK 



and (using (3.20) with k > and r > 0) 

o = 0, r 



K . 



(4.13) 



(4.14) 



Inserting (4.14) in (4.1 )-(4.4) and localizing the fields in the same edge (i.e. setting i = j) 
one finds that the correlation functions simplify to 



C n (vt 12 + x 12 , i; 0, fi) = (ipiih, xi, i)ipi(t 2 , x 2 , 



2tt V 2i 



3 i(utl2+X12)Mi 



sinh 



(vti2 + X12) - ie 



(4.15) 
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C 2 2(vti2 ~ Xi2,i] P,(j) = (ll>2{tl,X 1 ,i)lf)2{t2,X2,i))l3,n 

1 



hi:- 



v 7 k=i 



(vtl2—X12)S ik fJ,k 



sinh 



fa 



(vt 12 - x 12 ) - ie 



(4.16) 



The condition g + = g_ and the left-right asymmetry of the NESS construction in section 
2 imply that only left moving (incoming) excitations contribute to Cu, which therefore 
coincides with the equilibrium correlator [19]. All the non-equilibrium features are cap- 
tured by C22, which involves only right moving (outgoing) excitations. In fact, in spite 
of being localized in the edge E{ of the graph, (4.16) depends on the temperatures and 
chemical potentials of the all n edges. 



It is instructive for this reason to derive and compare the Fourier transforms of (4.15 



4.16). We will show first that they can be expressed in terms of the finite temperature 



TL any on distribution discovered in [19] . Consider in fact 

/oo /»oo 
dt / dxe-' l{Evt+px) C u (vt + x,i;(3,fx) . (4.17) 
-00 J —00 



Plugging (4.15) in (4.17) one gets 



7T 



C u (E,p,i;P,[i) = — 5(E-p) d(p - ft, fit k) 
v 



(4.18) 



where the 5-function fixes the dispersion relation and d is the equilibrium anyon momen- 
tum distribution flUl 



d(p,/3;K) 



pi-* e -()p/2 

pl-* e -pp/2 

2tvT(k) 



K 1 n K 1 n 



1 



2vr 



1 



2tt 



K > 



(4.19) 



B and V being the beta and gamma functions (Euler's integrals of first and second kind 
respectively). Notice that for k 7^ 1 the distribution (4.19) depends on both (3 and p and 
not only on the dimensionless combination pp. Eq. (4.19) defines a smooth function of 
peR, which satisfies 



lim d(j>, P; k) - 

p— >o 

and has the following asymptotic behavior: 



P 1 - k Y 2 {k/2) 
2tiT{k) 



lim d(p, P\ k) 

p— >oo 



lim d(p, P; k) 

p— >— 00 



0, V k > 0, 



00 



< K < 1 
K — 1 , 
K > 1 . 



(4.20) 

(4.21) 
(4.22) 
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For positive integer k (i.e. for fermions and bosons) the distribution (4.19) simplifies to 



(1+e^P) [2(n-l)}\ Hj=l W 



2n 



2n - 1 



(4.23) 



(ePp— 1) (2n— 



k = 2n . 



where n = 1,2,... and the familiar Fermi and Bose distributions appear as prefactors. 
The first two fermion and boson distributions are 



d(p,#l) 



(1 + e?P) ' 
1 p 
(e^ — 1) 27r 



d(p,P;3) 
d(p,/3;4) = 



(vr 2 +p 2 /3 2 ) 



;i + e^) 8tt 2 /3 2 
1 p(4vr 2 + p 2 /3 2 ) 
- 1) 48vr 3 /3 2 



(4.24) 
(4.25) 



As expected, in the fermion point k = 1 of the TL liquid one gets the familiar Fermi 
distribution. In spite of the fact that the remaining boson and fermion points (k = 2, 3, ...) 
have been established in [13 120] more then a decade ago, to our knowledge their physical 



meaning and potential applications of (4.23) have not been fully explored 




Figure 3: The distribution d at fixed temperature j5 = 1 for k = 1/2 (dotted blue line) 
and ft = 3/2 (dashed black line) compared to the Fermi distribution k = 1 (continuous 
red curve). 



we 



In order to give an idea about the anyon distributions in the interval < k < 1 
show some of them in Fig|3j where the standard Fermi distribution (continuous red curve) 
is given for comparison. Figji] displays the behavior of the anyon distribution (4.19) for 
fixed k — 1/4 and different temperatures. For < k < 1 and with decreasing of the 
temperature T ~ 1//3 one observes the formation of a sharp peak at p = (in agreement 



with eq. (4.20)), which signals a condensation-like phenomenon 

. of (4.16), it is useful to cor 

0), the system being driven away from equilibrium 



Concerning the Fourier transform of (4.16), it is useful to consider first the case when 
all the temperatures are equal (/% 
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2.0 
1.5 




Figure 4: The distribution d at fixed k = 1/4 for different temperatures (3 = 0.2 (continuos 
red line), (3 = 0.4 (dashed black line) and (3 = 0.8 (dotted blue line). 



only by the voltages V{. In this case 

C 22 (vt 12 ~ Xi2,i;/3,H) = (^2(*i,Zi,z)V>2(*2,Z2,«))/8,A« 

1 



-(- 

2tt V 2i 



3 i(vtl2-SBl2) Efc = i §ifcMfc 







sinh 



|(fti2 - x i2 ) - ie 



(4.26) 



and therefore 



t~ it 

_ m-/v 

C 22 (E,p,i : P,n) = — 6(E-p) dip- } S ik iJ, k ,/3;K) 



(4.27) 



fe=i 



One has still the equilibrium distribution, with the energy shifted by a linear combination 
of the chemical potentials Hk, whose coefficients are the S-matrix elements. 

Finally, in the coordinate space the general expression (4.16 ) is a product of Cn-factors 
with different temperatures and chemical potentials. One gets therefore in momentum 
space the nested convolution formula 



7T 



C 22 (E,p,i:(3,v) = — 6(E-p) 



~2tt 



2tt 



dk 



n-l 



2tt 



X 



d{h - Snfj-i, (3i, nS 2 a )d(k 2 - S i2 ft 2 - k u (3 2 ; k§? 2 ) ■ ■ ■ dip - § in /i n - fc n _i, /3 n ; «S? 



(4.28) 



Being a convolution of distributions, (4.28) is also a well defined distribution. The NESS 
has therefore a remarkable property: the associated non- equilibrium distribution is 
simply a convolution of equilibrium distributions with different temperatures and chemical 
potentials. 

Since the general form of (4.28) is quite complicated, it is instructive to consider below 
case n = 2 a 

D 2 (p; Pi,fo;K, 



the case n = 2 and fJ,i — pL 2 — 0, focusing on 

w dk 



2tt 



d(k, (3±; «; cos 8) d(p — k, /3 2 ; k sin 



0g[O,tt). (4.29) 



1* 



Using the x-space representation, at equal temperatures one finds the relation 

D 2 (p;f3 J f3;K,6) = d(k, P; k) , V 9 G [0, vr) . (4.30) 

For 0i 7^ 02 the convolution D 2 defines a new distribution. Since we were not able to 
determine its explicit analytic form, we give some plots which are obtained numerically. 
The plots in Figs. [5] and [6] illustrate the behavior of D 2 for different values of k and 9. 
We see that even for 0i ^ (3 2 the distribution D 2 is similar to d with a kind of "effective" 
temperature and statistical parameter depending on /?i j2 , K and 9. 



2.0 




Figure 5: The distribution D 2 at fixed K = 1/2 and 9 = tt/A for different temperatures 
(/5i, /3 2 ) = (1/2, 1) (continuos red line), (/5i, /3 2 ) = (1,2) (dashed black line) and ((3i, (3 2 ) = 
(2,4) (dotted blue line). 




Figure 6: The distribution D 2 at fixed 0\ = 1, (3 2 = 6 and 9 = n/Q for statistical 
parameters k = 1/3 (continuos red line), k, = 2/3 (dashed black line) and k — 4/3 (dotted 
blue line). 



Summarizing, we derived above the two-point TL anyon correlation functions away 
from equilibrium. The results ( 4.1|4.5 ) are expressed as products of S-dependent pow- 
ers of equilibrium correlators at different temperatures. In agreement with this fact the 
momentum space anyon NESS distribution is the convolution of equilibrium anyon distri- 
butions (4.28). The equilibrium limit satisfies the KMS conditions. At criticality one is 
dealing with a c = 1 conformal field theory, whose anomalous dimensions (4.12) depend 
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not only on the coupling constants g±, but also on the scattering matrix S. The above 
technique allows to compute higher anyon correlation functions as well, but in order to 
investigate the transport properties of the system, we concentrate below on the electric 
and energy current correlators away from equilibrium. 

4.2 Charge and heat transport 

The charge transport in the NESS is described by 



(dtQ+)p,p 



2< ; 



2ttv f (1 



J2( 5 i3 ~ > 



(4.31) 
(4.32) 



which follow by substituting (3.1313.28) in (3.34). Eq. (4.31 ) describes the external charge 



flow in the junction: it is constant in time and is incoming for (dtQ+)p tft > and outgoing 
for (d t Q+)p >ft < 0. Eq. (4.32) determines instead the value of the currents along the leads. 
The charge balance 



(4.33) 



is satisfied and represents an useful check. If S G O v , the electric charge is conserved 
(d t Q+)p )fi = and the fc^-dependence in (4.32) drops out, leading to 



(3 + (t,x,i))p u 



(4.34) 



3=1 



The current (4.34) satisfies the Kirchhoff rule (3.30) and vanishes at equilibrium (Vj = V 
for all i) as it should be. The dependence on the statistical parameter k is explicit and 



in front of the sum in (4.34) is positive 



deserves a comment. In the physical domain V, defined by (3.24), the overall coefficient 

TTKVp + 2g_ 



G(g-,K) 



2hvfC+ 2w 2 k 2 vf 



> 0. 



(4.35) 



For g_ > the coefficient G decreases monotonically with k > 0. For g_ < one has 
that k > —2g-/irvF in the physical domain T>. In this case the coefficient G increases 
in the interval —2g_/7TVF < K < —^g^/irvF, reaching the maximal value — v F /16g~ and 
decreases for k > —Ag_/ixvf. This behavior is illustrated in Fig. [7j 

The current (4.34) is proportional to the applied external voltages. Non-linear effects 



are absent in the critical regime under consideration, which implies the conductance tensor 

(4.36) 



GG7_,K)X;(Sii-a 



13 J 
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Figure 7: Behavior of G at = 1 (dashed black line) and = — 1 (continuous red line). 



We see that the NESS approach, adopted in this paper, confirms the result for G, obtained 
previously for k = 1 by different methods, including renormalization group techniques 
|H1 09] , linear response theory [26l 127] and conformal field theory [50| |5T] . The novelty 
in (4.36) is the explicit dependence on the statistical parameter k, shown in Fig. [7| 
A similar computation gives the energy (heat) flow 



(9 xt (t,x,i))^ 



to 



2 n I 



(4.37) 



which satisfies the Kirchhoff rule (2.26) for all § G 0(n). For § G O v the expression (4.37) 
takes the form 



n / n 

v? - 2k F J2&i - - 



Eft 

j=i 



«3 2 1 — 

12 ^/? 2 ' 



(4.38) 



The heat flow depends therefore not only on hp and the voltages V$, but also on the 
temperatures 

As before, we consider for illustration the case of n = 2 wires. Inserting the scattering 



matrices (2.31) in (4.32) one has 



Therefore, 



Gvp[(fii + /i 2 )(cos# — 1) — — // 2 ) sin^] , det§ = 1 , 
Gvp[(fii + /i 2 )(sin 9 — 1) + (/i! — // 2 ) cos 0] , det § = — 1 . 



(4.39) 



# = =>j+(t,x,l)=j+{t,x,2) = 0, 
= tt/2 x, 1) = -j+(t, x, 2) = G(ni - /i 2 ) 



det § = 1 , 
det§ = -1 , 
(4.40) 
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corresponding respectively to full reflection (disconnected edges) and complete transmis- 
sion in the junction. Finally, 



(9 xt (t,x,l))/3^ = -{6 xt {t,x,2))^ 



8ttC 



[(Hi — /x 2 .) sin 2 9 — /ii/i2 sin 29] 



irv 2 sin 2 9/1 



12 



2 ; 



(4.41) 



for both families in (2.31) 



4.3 Quantum noise 

In this section we derive the noise power in the TL junction in Fig. [TJ For this purpose 
we need [52] the two-point connected current- current correlator 



v conn 



(j + (t 1: x h i)j+(t 2 , x 2 , 
(j+(h, xx,i)j + (t 2 , x 2 ,j)) p>lt - (j+(h, xi,i))fr lt (j + (t 2 ,x 2 ,j))fo 

After some algebra one finds 

( ]+ {t 1 ,x l ,i)j + {t 2 ,x 2 ,j)) c ^ = 
V \ I 1 „ A„ 1 



(4.42) 



2v f(+J i/3 2 sinh 2 f.{vt X2 + x 12 ) - i: 

1 



7T^ + E S 



i = i [3? sinh |-(wti2 - x 12 ) - \e 
1 



/3| sinh 2 jfr(i;ii2 - £12) - ie 



/3 2 sinh 2 



^(u*i2 + aria) - ie 



J*. 



. (4.43) 



One easily verifies that the equilibrium limit (/% — >■ /3 for all i) of ( |4.43 ) satisfies the KMS 
condition 



(4-44) 

where the KMS automorphism £> acts on j + as follows, 

[Qsj+\(t,x,i) =j+(t + s/v,x,i). (4.45) 



The explicit expression (4.43) contains fundamental physical information about the 
NESS. First of all, since 



(j + (ti, x 1 ,i 1 )j+(t 2 , x 2 , i 2 ))™™ ^ (j+(-tx,Xi, ix)j + (-t 2 , x 2 , i 2 ))| 



conn 



(4.46) 



the NESS breaks down time reversal invariance, even if the junction interaction preserves 
it, i.e. if § = §* [47]. Nevertheless, time translation invariance is preserved, which allows 
one to use the conventional definition [52J of noise power 



/CO 
dte iuJt (j x (t,xx,i)j x (0,x 2 ,j))p 
-co 



conn 



(4.47) 



'The bar indicates complex conjugation. 
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Eq. (4.47) defines a complex matrix whose entries can be expressed [53j in terms of the 
hypergeometric function 2 F\, namely 



2v F (, 



{[F-(u,fr 



X12) 



F + {u,P u x 12 )]p- 1 5 lj + 



^2 ^ilPf 1 [ F - ( w > Ph -E12) - F+(w, Pi, -X12)] - 



1=1 



Si^lF-faP^-xu) -F+faP^-xu)] - [F_(u,p h x 12 ) -F+faPi^nftp^} 



with 



F±(u,p,x) 



e ±27Tx/l3 

iujp ± 2tiv 



(4.48) 



(4.49) 



From (4.48, 4.49) one can deduce the zero-frequency limit (zero-frequency noise power) 

P tj {P) = lim P^p; Xl ,x 2 ; u) . (4.50) 

Using 

(4.51; 



lim [F_(w; /3, x) - F + (u; p, x)} = 

w->0+ ZTTV 



one gets 



^f 1 U% - S^" 1 - K% + ^Pr%) 



(4.52) 



where G, defined by (4.35), captures the dependence (see FigQ of the noise on the 
statistical parameter k. As expected, Pij(fi) turns out to be a £i )2 -independent real 
symmetric matrix. If the electric charge is conserved (S G O v ), the noise power (4.52) 
satisfies in addition the Kirchhoff rule 



(4.53) 



i=i 



The expression (4.52) admits the typical Johnson- Nyquist (3" 1 behavior and shows the 
non-trivial interplay between the different temperatures and the scattering matrix. For 
example, in the two-terminal case with § = § + one finds 

G ( Ti(l -cos#) 2 + T 2 sin 2 # (T\ - T 2 )(l - cos#) sin# ^ 



" v F k B V (Ti-r 2 )(l-cos0)sin0 T 2 (l -cos^ + TiSin 2 ^ 

where T = (fc^/?) -1 is the absolute temperature and kg is the Boltzmann constant. The 
eigenvalues of P + 

1C 

pf = —r-i 1 -cos6)Ti > 0, i = l,2 (4.55) 



vpk 



B 



are nonnegative in agreement with the positivity of the two-point function (4.42). Anal- 
ogous result holds for P~ corresponding to S>~. 
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5 Outlook and conclusions 



In this paper we constructed and investigated an exactly solvable model of a non-equilibrium 
Luttinger junction. The basic points of our approach are: 

(i) a scale invariant point-like interaction, which is described by a scattering matrix § 

and drives the system away from equilibrium; 

(ii) a representation generated by a NESS fyg,^, which encodes the point-like interaction 

in the chiral fields (p^z', 

(iii) an exact operator solution of the TL model (in terms of tpi t z) on a star graph with 
the current splitting boundary condition in the vertex; 

(iv) an extension of the conventional fermion Luttinger liquid to anyon statistics. 

Combining these ingredients, we derived the basic correlation functions in the state Qp^- 
The essential characteristic features of these functions are: 

(a) the non-equilibrium two-point anyon correlations are products of §-dependent powers 

of equilibrium correlations at the temperatures and chemical potentials of the heat 
baths, connected to the leads; 

(b) accordingly, the corresponding momentum space distribution is the convolution of 

equilibrium anyon distributions at different temperatures and chemical potentials; 

(c) the Fourier transform of the leading terms in the large distance expansion of the anyon 

correlations gives Cauchy-Lorentz distributions, which after convolution reproduce 
themselves with appropriate width and median; 

(d) in the critical limit one has a c = 1 conformal field theory with S-dependent anoma- 

lous dimensions, which are explicitly derived; 

(e) the expected breakdown of time reversal invariance is manifest in the current-current 

correlator. 

We investigated in detail the energy and charge transport in the junction for all values 
of the statistical parameter. The energy is conserved for § e 0{n), which covers both 
possibilities of a junction without and with electric charge dissipation. In the latter case we 
determined the exact expression for the charge flow leaving or entering the junction. The 
connected current- current correlation is a linear combination of hypergeometric functions. 
The associated zero-frequency noise power depends linearly on the temperatures. 

Our investigation above has been focused essentially on the critical properties of anyon 
Luttinger liquids away from equilibrium. It will be interesting to study the noncritical 
aspects as well. The generalization of the results of this paper beyond the Luttinger 
liquid paradigm, when the nonlinearity of the dispersion relation becomes essential, is 
also a challenging open problem. 
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A Large space separation asymptotics 



The behavior of the correlators ( |4.15|4.16 ) at large large space separation |x 12 | 3> Pi is 
encoded in 



C u (x,i; (3,fi) 
C 22 (x,i;(3,jj) 



r / ii i ) I L(x;^^,-§ ifc /i fc ) + 



1 \ -i— r / Zix 



n 



(A.l) 
(A.2) 



where x 3> Pi > 0, the dots stand for sub- leading contributions and 

1 



L(x; 7 ,/i) = -^e" 1 ^-^, x>0. 

Z7T 



The Fourier transform 



dxe ipx L(x;7,/i) 



1 



7 



VT7 2 + (p - /i) S 



(A.3) 



(A.4) 



is the familiar Cauchy-Lorentz distribution^} where 7 is the half width at half maximum 
and fj, is the statistical median. Using that the class of Cauchy-Lorentz distributions is 
closed under convolution, one finds 



Cu(p,z;/3,/i) 
<? 22 (p,z;/3,/i) 



dxe ipx C n (x,i; p,n) ~ L ( p; — , 



7T/t 



(A.5) 



n s 2 . 



dx e iprr C 2 2(a;, i; P, l*) ~ L ( p; ttk y] - y~] Sy/ij ] + 



(A.6) 



Summarizing, the Fourier transform of the leading term in the long distance expansion 



of both (4.15) and (4.16) is a Cauchy-Lorentz distribution. Notice that the width and the 



median of (A.6) depend on the temperatures and chemical potentials of all heat baths, as 
well as on S. 



7 Known also as non-relativistic Breit-Wigner distribution. 
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